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Abstract The Dirac equation is solved in the rotating Bertotti-Robinson 
spacetime. The set of equations representing the Dirac equation in the 
Newman- Penrose formalism is decoupled into an axial and angular part. 
The axial equation, which is independent of mass, is solved exactly in terms 
of hypergeometric functions. The angular equation is considered both for 
massless (neutrino) and massive spin-^ particles. For the neutrinos, it is 
shown that the angular equation admits an exact solution in terms of the 
confluent Hcun equation. In the existence of mass, the angular equation 
does not allow an analytical solution, however, it is expressible as a set of 
first order differential equations apt for numerical study. 
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1 Introduction 

The Rotating Bertotti- Robinson spacetime (RBR), which was discovered 
long time ago by Carter [Tj, is an Einstein-Maxwell solution representing a 
rotating electromagnetic field. This solution remained unnoticed until the 
study of Al-Badawi and Halilsoy [2], who rediscovered it by applying a 
coordinate transformation to the cross-polarized BcU-Szekeres solution of 
colliding electromagnetic waves [3l. We can consider the RBR solution as 
an extended version of the Bertotti- Robinson (BR) solution due to the fact 
that the RBR solution contains one more degree of freedom to be interpreted 
as rotation. Adding rotation to the BR creates gravitational curvature, dis- 
torts isotropy and modifies geodesies significantly. For this reason the RBR 
solution assumes more complicated topology compared with the BR solu- 
tion. The RBR spacetime has the topology of AdS^ x and underlying 
group structure of SL{2,R) x C/(l). Nowadays, the spacetimes with AdS 
structure are quite popular because of their connection with string theory, 
higher dimensions and the brane worlds. The RBR solution can also be in- 
terpreted as the "throat" connecting two rotating black holes with charges. 
This is due to the fact that the BR solution is considered as the "throat" 
connecting two asymptotically flat Reissner-Nordstrom regions [4] . 

The first study on the Dirac equation in the BR spacetime without 
charge coupling was considered long time ago |5j. During the last decade 
the studies on spin-i particles in the BR spacetime gained momentum. For 
example, Silva-Ortigoza |6] showed how the Dirac equation could be sepa- 
rated when the background is the BR spacetime with cosmological constant. 
Later on, a more detailed study on the problem of the Dirac equation in 
the BR spacetime was worked out as well [TJ. Here, we extend this recent 
work by looking for the answer to the following question : "How does a 
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test Dirac particle behave in a rotating spacetime filled with electromag- 
netic field i.e., in the RBR spacetime?". We shall ignore the back-reaction 
effect of the spin-i particle on the spacetime, by the same token as done 
in [7] . The RBR solution represents one of the type-D spacetimes and as it 
could be followed from the literature the studies on spin-i particles in type- 
D spacetimes have always attracted attention 8 (and references therein). 
More recently, the problem of Dirac equation in the near horizon geometry 
of an extreme Kerr black hole (Kerr throat) has been studied by [9]. In many 
aspects, that spacetime of the Kerr throat shares common features with the 
RBR spacetime. The main difference between them is that the Kerr throat 
is a vacuum solution, while the RBR is not. On the other hand, they are 
both regular solutions. 

In this paper, our aim is to solve the Dirac equation in the RBR space- 
time. To this end in order to separate the Dirac equation we employ the 
well-known method suggested by Chandrasekhar [8j . We separate the Dirac 
equation into the axial (function of z only) and the angular parts (function 
of d only) in such a way that the resulting axial equation remains inde- 
pendent of mass. This advantage leads us to obtain an exact solution for 
the axial equations in terms of hypergeometric functions. The angular part 
turns out to be more complicated than the axial part. This is due to the 
fact that the metric functions are dependent on the variable 0, and also the 
angular equations contain the mass term. For the angular part two separate 
cases, massive and massless (neutrinos) cases are discussed. The angular 
equations of the massive case, we are able to reduce the equations to a set 
of linear first order differential equations, which can be utilized numerically. 
However the massless particle (neutrino) equations are solved exactly by 
reducing the equations to the confluent Hcun equations. 
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Recall that, the confluent Heun differential equations are less known 
than the hypergeometric family in literature. The modern mathematical de- 
velopment shows that many physical problems are solved exactly by Heun 
functions. |10 p i H [T2]. For example, problems involving atomic physics with 
certain potentials I13J which combine different inverse power or combin- 
ing the quadratic potential with inverse powers of two, four, and six etc.. 
Problems in solid state physics, like dislocation movement in crystalline 
materials and quantum diffusion of kinks along dislocations are also solved 
in terms of Heun function [14 . Problems in general relativity, Fiziev [15| 
gave an exact solution of the Regge- Wheeler equation in terms of the Heun 
functions and apply them for study of different boundary problems. More 
recently, Birkandan and Hortacsu |16j have given examples in which the 
Heun functions admit the solution of the wave equation encountered in 
general relativity. They have related the solutions of the Dirac equation, 
when the background is the Nutku's helicoid spacetime in five dimensions 
to the double confluent Heun function. Nowadays, the modern computer 
packages have started to involve the Heun functions in their algorithms, as 
for instance it can be seen in the 10*'* and higher versions of the famous 
computer package. Maple. 

The paper is organized as follows: In section II, a brief review of the RBR 
solution is given. Next, we present the basic Dirac equations and separate 
them in the spacetime of RBR. In section HI, we present the exact solution 
of the axial equation. The angular equation with both massless and massive 
cases are discussed in sections IV and V, respectively. Finally, in section VI, 
we draw our conclusions. 
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2 Rotating Bertotti-Robinson spacetime and separation of Dirac 
equation on it 



The metric describing a rotating electromagnetic field, RBR solution, writ 
ten in spherical coordinates, is given by [2] 

m 



dp - dr^ - r^de^ 



sin^ I 



F2 



-dt 



(1) 



where the function F{9) and the constant q are 

F{e) = 1 + (1 + cos^ e) , 

q = 2a\/l + a^. 



(2) 



in which a is the rotation parameter. It is readily seen that for a = 
metric (1) reduces to the BR metric [7] (and references cited therein). 

We make the choice of the following null tetrad basis 1-forms (Z, n, m, m) 
of the Newman-Penrose (NP) formalism [17J in terms of the RBR geometry 
that satisfies the orthogonality conditions, l.n = —m.rn = 1. We note that 
throughout the paper, a bar over a quantity denotes complex conjugation. 
We can write the covariant 1-forms as 

^/2l ^ —^/Fidt- dr) , 
2r 



\/2ri = -VF{dt + dr) , 

r- ,„ sin 9 /2a / -~ 

V2m = iVFde + —== f — VI + a^dt 



^ . V . _ , . (3) 

VF \r ^ ' ^ ^ 

We obtain the non-zero 'F2 and <Pii, which are known as Weyl and 
Maxwell scalars, respectively, as 

1'2 = ^ 



F3 



(1 + a^) cos 29 + cos^ 9 vT+o^ (1 + + sin^ 9) cos t 
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The singularity free and the type-D characters of the metric can be 
easily deduced from ^2- It is obvious that for a = 0, this typc-D metric 
(1) yields a conformally flat spacctimc (i.e. the BR spacctimc) in which a 
uniform electromagnetic field, with = i, fills the entire space. As it 
can be seen from Eq. (4), rotation (a ^ 0) gives rise to anisotropy of the 
prevailing electromagnetic field. 

In order to study the Dirac equation in the RBR spacetime we prefer 
to work in a more convenient coordinate system, therefore by using the 
following transformation. 



t = tan 



2t 



metric (1) is transformed into 



= F{e) 



de^ 



{d(f> — qzdty 



(5) 



(6) 



(1+^2) f2 

We notice that the metric functions in Eq. (6) depend explicitly on 
the variable 6 as a result the angular part of the Dirac equation becomes 
important. The coordinates — oo <t<oo, — oo < 2: < oo covers the entire, 
singularity free spacetime. 

The covariant 1-forms of the metric (6) can be taken as. 



V2I = VF{y/lTz^dt ■ 



dz 



dz 

7f+7 
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V2m = y/Fde 



I Slllt 



-.{d(j) — qzdt), 



while their corresponding directional derivatives become 
dt 



V2D 



V2A 



FVl 
dt 



V2S- 



V2S- 









fF 








fF 


* + 


.Vf 




smb' 


dg 






I 

sinfc' 



qzd^ 



(7) 



(8) 



Using the above tetrad we determine the nonzero NP complex spin co- 
efficients 117' as, 



-1 



T = — TT = 



2V2FI 



[a^sin(26') + igsin6')] , 



e = 7 



a = -/3 



[2cot6'(l + 2a^) + iqsinl 



The Dirac equations in the NP formalism are given by [8^ as 

{D + e- p)Fi + (S + n - a) F2 ^ ifipGi, 

(5 + /3 - r) Fi + (Z\ + /i - 7) F2 - iMpGs, 

{D + e-p)G2- (S + Tf -a)Gi = i/i^Fj, 

(Z\ + /I-7) Gi- {5 + P-f)G2 =ifipFi, 

where /i* = V^Pp is the mass of the Dirac particle. 

The form of the Dirac equation suggests that we assume [8], 

Fi = (61) 6^^=*+"^^ 



(9) 



(10) 
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F2 



/a {z) {6) e'(*^*+™'^) 



51 {z)A2 {0) e^('=*+'"<^), 



G2 



52 {z)Ai {9) e*(*=*+™'^). 



(11) 



Here, we consider the corresponding Compton wave of the Dirac parti- 
cle as in the form of / (2;) A (9) e'(*'*+™'^) , where k is the frequency of the 
incoming wave and m is the azimuthal quantum number of the wave. The 
temporal and azimuthal dependencies are chosen same but axial and angular 
dependencies are chosen different for different spinors. 

Substituting the appropriate spin coefficients (9) and the spinors (11) 
into the Dirac equation (10), we obtain 





Zg2 ^ £^A2 . J2A3 r- 

1 -. = ILl -T-vF. 



91 Ai gi A4 



Z91 £Ai 




(12) 



91 M 



where the axial and the angular operators, respectively, are 



Z = Vl + z^d, + 



1 



[z + 2i{k + mqz)] , 



and 



Z = Vl + z^d, + ^ i [z - 2i{k + mqz)] , 

2vl + 2; 
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^ + „ cot 9 0? sin 2Q mF q sin ( 
=de + — — — - i- 



4F sin 6* 4F 



cot sin 2Q mF . g sin ( 
£ = ^0^ -— h + 1- 



4F sine 4F ' 

_ cote sin 20 mF asm.Q 
£^ = de^-^ ^+^^- 

It is obvious that L and -L^ are purely angular operators and L = £, 

In order to separate the Dirac equation (12) into axial and angular parts 
we choose /i = 52, /2 = 51, ^2 = ^1 and = As and introduce a real 
separation constant A as 

Zg2 = -Xgi, (14) 



Zgi = -A52. (15) 



and 



LA3+ifi*AiVF = -XAi, (16) 
Vf = XA3. (17) 

3 Solution of the axial equation 

The structure of the axial Eqs. (14, 15) admits gi = Thus it is enough 
to decouple the axial equations for g2, namely 

I (Z52) = X''g2, (18) 
The explicit form of Eq. (18) can be obtained as 
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(1 + z^)g'^ {z) + 2zg'2 (z) + 



2(1 + ^2) 



1 + :^^ + 



2{k + mqzf - 2i{kz - mq) - 2A^(1 + z"^)] {z) = 0, (19) 

(Throughout the paper, a prime denotes the derivative with respect to 
its argument.) 

Let us introduce a new variable y such that z = i{l — 2y), therefore Eq. 
(19) becomes, 

y{l-y)g'i{y) + {l-2y)g'^ [y)- 



42/(1 



^1- 

-y) 14 



mV(l-2y)2 + 2/(l-y)+ 



k{k + l- 2y) - 4A\(1 -y)+ imq [1 + 2fc(l - 2y)]} g^ {y) = 0, (20) 

The exact solution of the axial part is found in terms of the Gauss 
hypergeometric functions as, 

52(y) =Ciy"(y-l)'5 2Fi{^+k-'y,^ + k + 'y,^+k + imq;y)+ 

C2y^"{y- 1)^ 2Fi{-{'y + imq),j - imq,^ -k-imq;y), (21) 

where the parameters are 

a=i(t + i + .™,). 

/3= i(fc- i-img). 
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7 = -ra^q^. (22) 
and Ci , C2 are complex constants. 

4 Reduction of the angular equation to Heun equation: The 
massless case 

The aim of this section is to show that the angular Eqs. (16, 17) for the 
massless Dirac particles (such as neutrinos) can be decoupled to a confluent 
Heun differential equation. 

For /i* = 0, Eqs. (16, 17) can be explicitly written as 

{Q) + {K-M) Ai {9) = XA3 (0) , (23) 
A'^{e) + {K + M)A3{e) = -XA,{e), (24) 

where 



cot (9 iiin20 qiiin0 , , 
(2^) 



M=^, (26) 
smt* 

Introducing the scalings 

Ai {0) = i?i {0) e- n^-^)d0^ (27) 
As {0) = Ha {6) e- (28) 

we get 

H[{0) = XH3{0)e-^^^'''^, (29) 
ff^ {0) = -XHi {0) e'^f^'i^, (30) 
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Decoupling Eqs. (29, 30) in Eq. (29) for Hi (0) , we obtain 

H'^ (9) + 2MH[ {6) + X^Hi (9) = 0, (31) 

Similarly, if we decouple the axial equations for H3 (9), the resulting 
equation turns out to be 



H'^ [9) - 2MH'^ {9) + A^i/g {9) - 0, 



(32) 



Introducing a new variable 9 = cos ^ (1 — 2z) , Eqs. (31, 32) are cast 
into the general confluent form of Heun equation, namely 



T i + m + 2a^m ^ — (ni + 2a^m) 
'Aa^m + h - 



z - 1 



„ i — (m + 2a?m) ^ + m + 2a?m 
Aa^m + ^ ^ + ^ 



z-1 



H[iz)--^^miz)^0, 
(33) 

H:,{zy^^H,iz) = o. 



z{z^l) 



(34) 



Recall that the general confluent form of Heun equation [T^] is given as 
follows 



H"{z)- 



A 



B C 



z z — 1 



n'{z)^^^B{z)=Q^ (35) 



The 



After matching Eqs. (33, 34) with Eq. (35), one can get the following 
correspondences . 
a) For Eq. (33), 



A = -4a^rn, B = ^+m+2a^m, C = i-(TO+2a^m), £» = and /i = A^. 



(36) 



b) For Eq. (34), 



A = 4a^TO, S = i-(m+2a^m), C = i+TO+2a^TO, D = and /i A^. 

(37) 
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Determining when the solutions of a confluent Heun equation are ex- 
pressible in terms of more familiar functions would obviously be useful. 
Expansion of solutions to the confluent Heun equation in terms of hyper- 
geometric and confluent hypergeometric functions are studied in detail by 
[l2] . In [12], it is also shown that the confluent Heun's functions can be 
normalized to constitute a group of orthogonal complete functions. Here, 
as an example, we follow the intermediate steps in [12] (page 102) in order 
to express the solutions of Eq. (35) with Z? = in terms of the hypergeo- 
metric functions. The transformation between the confluent Heun function 
and the hypergeometric function is given with the Floquet expansion [12[ 
[T4] . Namely, 

OO 

H,{z)= 5n 2^1(^1,^2; 5; z), (38) 

n— — 00 

where 

^1 = -71 - i^j and ^2 = n + I'j + C + B - 1, (39) 

and i^j are known as the Floquet exponents. The coefficients g„ satisfy 
a three-term recurrence relation 

Angn-l + QnQn + = 0, (40) 

where 

An = Abn-l,n, 
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•Pn = Abn+l,n, (41) 

and the coefficients, 6„_i_„, 5„_„ and 6„+i,„ expressed in terms of the 
parameters Cii?2i given explicitly by 42j .(The only difference between 
our notation and the Ronveaux's notation is -B = 7 ). 

Finally, it should be noted that one can see a brief analysis of the conflu- 
ent Heun equation as well as with its power series solution and polynomial 
solution in [18] . 



5 Reduction of the angular equation into a set of linear first 
order differential equations: The case with mass 

In this section, we shall reduce the angular Eqs. (16, 17) into a set of linear 
set of first order differential equations for the case of Dirac particle with 
mass. To this end, let us make the following substitutions into Eqs. (16, 
17), 



cot B a-' 5in(2t)) 

A^{e)^[A^{e)^iBa{e)]e '\ ' (42) 

A3 (0) = [Afo (0) + *A^o m e n ^ ) ^ (43) 

where Aq (6) , Bq (9) , Mq (0) and Nq (6) are real functions. After sepa- 
rating the real and the imaginary parts, one obtains a set of first order 
differential equations 

A'o (0) - ^Ao (0) + (0) = \Mo (9) - f^^VPN^ (6) , 

sm 4:1' 

B'o id) - —,B, {9) - i^Ao {e) = XN, {9) - fi*VFMo (9) , 
sm 9 At 

M'o (0) + (9) + ^^No (9) = -XAo (9) - h*VfBo (9) , 

sm At 
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K W + {6) - {6) = -XBo {6) - „*^/FAo {9) . (44) 

Introducing a new variable x = cos 6, one can remove the trigonometric 
functions in the set of first order differential equations (45). But solving 
the entire system analytically does not seem possible. To our knowledge, in 
the literature such a system does not exist. Nevertheless, one can analyze 
the system via an appropriate numerical technique, which may need an 
advanced computational study. 

6 Conclusion 

In this paper, our target was not only to separate the Dirac equation for 
a test spin- 5 particle in the rotating electromagnetic spacetime (RBR), 
but to explore exact solutions as well. By this way, we wanted to make a 
contribution to the wave mechanical aspects of the Dirac particles in the 
RBR geometry. 

Due to the metric functions, which are only functions of the angular vari- 
able 6, the angular part of the Dirac equation in the RBR background is the 
harder part to be tackled compared with the axial part. Another advantage 
of the axial part is that the axial equations do not involve the mass term. 
These simplifications in the axial equations guided us to obtain the general 
solution of the axial part in terms of the hypergeometric functions. On the 
other hand, although we could not obtain the general analytic solution of 
the angular part, we succeeded to overcome the difficulties in the angular 
part in the massless case and obtained an exact solution of it in terms of 
the Heun polynomials. Inclusion of mass prevents us to obtain an analytic 
solution for the angular part. As an alternative way to the analytic solu- 
tion, in the last section we showed that the angular part could be written 
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as a set of first order differential equations, which are suitable for numerical 
investigations. 

Finally, the study of the charged Dirac particles in the RBR spacetime 
may reveal more information compared to the present case. This is going to 
be our next problem in the near future. 
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